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Grover's quantum (search) algorithm exploits principles of quantum information theory 
and computation to surpass the strong Church-Turing limit governing classical computers. 
The algorithm initializes a search field into superposed N (eigen)states to later execute 
nonclassical "subroutines" involving unitary phase shifts of measured states and to 
produce root-rate or quadratic gain in the algorithmic time (0(A/^/^)) needed to find some 
"target" solution m. Akin to this fast technological search algorithm, single eukaryotic 
cells, such as differentiated neurons, perform natural quadratic speed-up in the search 
for appropriate store-operated Ca^+ response regulation of, among other processes, 
protein and lipid biosynthesis, cell energetics, stress responses, cell fate and death, 
synaptic plasticity, and immunoprotection. Such speed-up in cellular decision making 
results from spatiotemporal dynamics of networked intracellular Ca^+-induced Ca^+ 
release and the search (or signaling) velocity of Ca^+ wave propagation. As chemical 
processes, such as the duration of Ca^+ mobilization, become rate-limiting over interstore 
distances, Ca^+ waves quadratically decrease interstore-travel time from slow saltatory 
to fast continuous gradients proportional to the square-root of the classical Ca2+ diffusion 
coefficient, D^^^, matching the computing efficiency of Grover's quantum algorithm. In 
this Hypothesis and Theory article, I elaborate on these traits using a fire-diffuse-fire 
model of store-operated cytosolic Ca^+ signaling valid for glutamatergic neurons. Salient 
model features corresponding to Grover's quantum algorithm are parameterized to meet 
requirements for the Oracle Hadamard transform and Grover's iteration. A neuronal 
version of Grover's quantum algorithm figures to benefit signal coincidence detection 
and integration, bidirectional synaptic plasticity, and other vital cell functions by rapidly 
selecting, ordering, and/or counting optional response regulation choices. 



Keywords: biotechnology, calcium-induced calcium reactions (CICRs), cellular decision making, classical and 
quantum computation, inositol 1,4,5-trisphosphate receptors (IP3RS), intracellular calcium, neuronal plasticity, 
quantum molecular networks and memory 



INTRODUCTION 

Modern analog, digital, and quantum descriptions of phyloge- 
neticaUy diverse cell functions (e.g., Monod and Jacob, 1961; 
McAdams and Shapiro, 1995; McAdams and Arkin, 2000) date 
to Twentieth-century revelations in computational and infor- 
mation sciences (e.g., Szilard, 1929; Turing, 1936; Shannon, 
1938, 1948a,b; Landauer, 1961; Feynman, 1982; Deutsch, 1985). 
Continued advances in systems biology, synthetic biology, and 
micro- and nanobiotechnology increasingly drive states-of- 
knowledge and -art in computational cell biology toward trends 
in logic gate, circuit, and algorithm designs (e.g., Ehrenfeucht 
et al, 2003; Amos, 2006; Baumgardner et al, 2009; Friedland 
et al, 2009; Adamatzky, 2010; Clark, 2010a,b,c,d, 2011, 2012b, 
2013a; Norris et al, 2011; Karafyllidis, 2012; Mehta and Schwab, 
2012; Daniel et al, 2013; Goiii-Moreno et al, 2013; Ji et al., 
2013), especially for "programmable" group and solitary cellular 
decisions mediated by genetic, epigenetic, and somatic regulatory 
networks. Unsurprisingly, given their preeminent status as com- 
putational units (cf Koch and Segev, 2000; Grillner, 2006), single 
neurons are stOl favored models for bioinspired smart technolo- 
gies (e.g., Liu et al., 2013). Yet, despite technological interests in 



neuronal information processing attributes, serious application 
of quantum computational approaches toward study of adaptive 
cybernetic-like neuron behavior and physiology remains disap- 
pointingly slow, except as it may broadly relate to more-or-less 
controversial debates over the statistical mechanics nature of con- 
sciousness, decision making, and other psychological states and 
functions of humans and animals (cf. Beck and Eccles, 1992; 
Hameroff, 1994, 2012; Tegmark, 2000; Schwartz et al, 2005; 
Khrennikov, 2009; Pothos and Busemeyer, 2013a,b). The strange 
properties of quantum mechanics, such as superposition, entan- 
glement, interference, and tunneling (Boxl), can be harnessed 
to enhance the information storage capacity, processing speed, 
and fault tolerance of man-made computational systems (cf 
Nielsen and Chuang, 2000). Over recent decades, quantum infor- 
mation theorists have steadily identified and adapted quantum 
computational constructs believed to outperform the classical 
universal Turing machine (Feynman, 1982; Deutsch, 1985), sup- 
ported by the strong Church-Turing thesis to be the upper limit 
for powerful computational devices obeying classical informa- 
tion theory and physicochemical laws. Actual physical quantum 
computers are only now moving out of proof-of-concept stage 
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Box 1 I Glossary of terminology. 



Bell Basis States or Einstein-Podolsky-Rosen Pairs 

Named respectively after John Bell and Albert Einstein, Nathan Rosen, and Boris Podolsky, these four orthonormal quantum states, span- 
ning the two-qubit state space, form an essential computational basis for many fundamentally useful applications of quantum information 
theory. The states or pairs are: \ f)oo = (I00> + |11))/2V2, = (|oi> + |10»/2V2, \^)^^ = (|00) - |11»/2V2, = (|01> - |10»/2V2. 
Quantum entanglement between the two qubits of a state permits secure cryptographic protocols, such as quantum teleportation, and 
compressed information encoding and storage, such as superdense coding and quantum memory. 

Church-Turing Limit 

Upper computational bound of efficiency for classical computers independently determined by Alonzo Church and Alan Turing. The 
Church-Turing limit emerges from the Church-Turing thesis or conjecture, which equates functions computable on a Turing machine with 
those computable by an algorithm. The strong limit is believed to subtend allowable complexity of computations performed by quantum 
computers, as captured in David Deutsch's rigorous conceptualization of universal quantum computers. 

Eigenstates, Eigenvalues, and Eigenspaces 

An eigenstate or eigenvector, \ v}, is a nonzero state in a state space operated on by a linear function L, so that L\v} =ci>\v} with com- 
plex eigenvalue co. Eigenstates and eigenvalues are ascertained from the characteristic function, c{"a) — diet\L - Xl\, where det is the 
determinant function. The eigenspace of a>, a subset of the state space on which L acts, is the set of eigenstates with the eigenvalue co. 

Entropic Uncertainty Principle 

Information or entropy principle first derived by Iwo Bialynicki-Birula and Jerzy Mycielski and by David Deutsch from Werner Heisenberg's 
uncertainty principle. The principle is formally expressed in the strong condition as HIQ) + H(R) > 2 log2[1 /f(Q, /?)], where H(Q) and H(R) 
are the Shannon entropies of respective spectrally decomposed measurements Q and R of quantum state \ f) with probability distributions 
p(q) and p(r) and maximum fidelity or inner product f[Q, R) = maxqj\{q\r)\ between eigenvectors \q) and \r). As with Heisenberg's version 
involving the standard deviation of observables, the entropic uncertainty principle places an upper bound on attainable knowledge about 
quantum systems. 

Hadamard Transformation 

An important step for quantum parallel computing, the Hadamard transformation applies the Hadamard gate n times to n input qubits 
to initialize the data register of a quantum information system into superposition, so that multiple values of an index integer x can be 
simultaneously analyzed by a single function f. Common notation for the Hadamard transform is H®". 
Hermitian Operator 

A linear operator, also known as a self-adjoint operator, imposed on a vector space V. For a Hermitian operator L, there is a unique linear 
operator acting on V, so that an adjoint or conjugate vector \v)^ exists for every \v) in V. That is, L = L\ \v) = \v)\ and LL^ = L^L when 
L is normal. 

Landauer's Principle 

Principle postulated by Rolf Landauerto define the relationship between energy and computation. Landauer improved earlier ideas of Le6 
Szil^rd, John von Neumann, and other theorists to concretize the minimum amount of energy/information consumed during irreversible 
operations. Although Landauer's principle applies generally to energy/information dissipated as heat from work, it is usually placed into 
the context of memory erasure. For example, for a biological or technological computer with exhausted finite memory capacity, it is 
necessary to erase information for further computations. Landauer's principle states at least k^Tlnl of energy, where kg is Boltzmann's 
constant and T is ambient temperature in degrees Kelvin, must be transferred to the environment for erasure of one bit of information. 
The corresponding entropy rendering of Landauer's principle may be written without variable T. 

Quantum Gates and Circuits 

Logic gates and circuits constructed to perform operations based on quantum mechanics and information theory. Popular unitary single- 



qubit gates in matrix form include, for instance, the quantum Identity / = 



1 0 
0 1 



0 -/ 

/ 0 



, and quantum Flip or Pauli Z = 



gates, the Hadamard H — 1/2 



1/2 



gate, the quantum NOT or Pauli X - 



gate, and the Phase S = 



1 0 
0 / 



0 1 

1 0 



Pauli Y = 



gate. These and 



other gates may be assembled into quantum circuits, such as multiple-qubit controlled-NOT, controlled-Phase, controlled-Swap or Fredkin, 
and Toffoli gates, also used to transform input qubits. The above single-qubit gates establish with other quantum gates a discrete subset 
of logical primitives (i.e., gates and/or circuits) capable of unitary transformation and of emulating any other transformation to approximate 
computational universality. 

Quantum Mechanical Properties 

Statistical wave-particle features of quantum mechanical systems, such as quantum superposition, entanglement, interference, and tunnel- 
ing, not observed for deterministic classical Newtonian physical or Shannon informational systems. Quantum superposition describes the 
linear combination or addition of state or vector solutions to Schrodinger's wave equation (or other quantum state equation permutations). 
When superposed states are indistinguishable, they are said to be entangled. Quantum interference is the disruption of state or vector 
(e.g., a wave or particle) spacetime trajectories. And quantum tunneling is the transition of one state to another without surmounting clas- 
sical energy barriers required for transformation in classical physics. Such probabilistic effects are useful for development of information 
technologies and additional purposes. 

Quantum Networks 

Technological and biological networks whose connectivity tends to obey either Bose-Einstein or Fermi-Dirac quantum statistics rather than 
classical Maxwell-Boltzmann statistics. The strength of each network node is described as a separate fitness or energy level and nodal 

(Continued) 
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Box 1 I Continued 

links take on the identity of particle states functioning under associative-like preferential attachment rules. In such cases, control parameter 
T (i.e., local absolute temperature), which dictates system behavior, is often replaced with a computational annealing parameter, such as 
space, time, or the "critical tunneling field strength." Quantum networks may display the network analogs of Bose-Einstein condensation 
and the Pauli exclusion principle depending on statistical parameters. In addition, rate of state transitions or computational decisions in 
a quantum network follows nonlinear first-order Arrhenius kinetics associated with quantum tunneling, also making it a computational or 
network analog of the physical phenomenon. These properties of quantum networks strongly compare with networks capable of certain 
associative forms of learning and memory, such as Hebbian-type learning dependent upon mutually weighted nodal or synaptic strengths. 

Schrodinger Wave Function 

A solution to Erwin Schrodinger's wave equation used to describe the statistical nature of eigenstates that exhibit wave-particle duality. 
The wave function gives the probable energy of a wave in a 4D spacetime interval. Each wave is associated with a respective wave 
number related to wavelength or an energy level known as a quantum. Because the wave equation is additive, the distribution of wave 
numbers or energy levels occurring over the spacetime interval may be combined into a mixed or superposed wave packet. Pure states 
are represented by a single wave number or quantum. 

Unitary Operator 

A mathematical operation performed on a state space that satisfies the identity operator by being normal, having a spectral decomposition 
(i.e., capable of being reduced to additive components), and preserving the inner product of two vectors, such as a unitary phase shift of 
180° from state or vector |0) to |1 ) or |1 ) to |0) residing on the unit circle or the unit 3D sphere — the Bloch sphere. 

Universal Turing IVIachine 

An idealized computational machine with unbounded memory belonging to a class of devices introduced by Alan Turing, the universal 
Turing machine is one of the most powerful classical computers conceived and serves as model to examine issues of computational 
complexity. All Turing machines are prototype modern programmable computers capable of executing algorithmic routines of different 
complexity. They consist of four essential components: (1) a microprocessor-like finite state control that coordinates computing action, (2) 
a program, (3) a memory tape, and (4) a read-write tape head that points to the tape location currently accessible for read-write operations. 
Unlike other Turing machines, the only variable component maintained by the universal Turing machine is the initial contents of the tape. 
Such a configuration is deceivably powerful and enables the universal Turing machine to emulate or simulate the processing of all other 
Turing machines, even more powerful ones. 



due to gradual progress in innovating suitable, if not optimal, 
device architectures, such as optical lattices, ion traps, nuclear 
magnetic resonators, quantum dots, and other technologies (e.g., 
Vandersypen et al, 2001; DiCarlo et al, 2009; Politi et al, 2009; 
Johnson et al, 2011). Nonetheless, algorithms built from quan- 
tum gates and circuits offer exciting practical, though often intu- 
itively difficult, ways for obtaining performance characteristics 
better than those exhibited by classical processors. For example, 
several general classes of quantum algorithms based on Shor's 
quantum Fourier transform (Shor, 1994), the Deutsch-Jozsa 
algorithm, and Grover's quantum algorithm (Grover, 1996) are 
known. Through a series of quantum logic gates, Shor's quantum 
Fourier transform enables a fast two-register eigenvalue phase- 
estimation procedure to be executed on eigenstates of a unitary 
operator put into quantum superposition. Phase-estimation sub- 
routines, in turn, serve as modules for other algorithms that 
exponentially decrease the number of operations required to solve 
important problems, such as related order-finding and hidden- 
subgroup problems, judged intractable with classical computers. 
Alternately, Grover's quantum algorithm, also termed Grover's 
quantum search algorithm, the quantum search algorithm, or 
the fast search algorithm, reaches only root-rate or quadratic 
operating improvements when compared to classical algorithms 
searching or counting elements of unstructured databases. This 
single-register algorithm initializes the search field containing tar- 
get solutions into a uniform superposition state via a quantum 
transformation. A quantum subroutine called the Grover's oper- 
ator or iteration then conditionally shifts or rotates the phase of 
certain computational bases until search solutions become found. 



Similar processing advantages arising from quantum mechan- 
ics may exist for natural computations performed by live bio- 
logical systems, particularly at the level of single cells and their 
subcellular components. Mounting evidence from decades of 
analytical and experimental research continues to oppose the 
conventional tenet that quantum mechanical phenomena exert, 
at best, trivial influences over bioprocesses (cf Davies, 2004). 
Criticisms still tend to concentrate on the capacity of biological 
systems to settle or cohere into a quantum regime long enough 
to accomplish quantum computation (cf Tegmark, 2000; Davies, 
2004; Reimers et al, 2009; Wolynes, 2009; Trevors and Masson, 
2010). However, issues regarding quantum decoherence, the col- 
lapse of the Schrodinger wave function into a single classical 
or macroscopic state due to thermodynamic processes involving 
a system and its environment, are less problematic for cellular 
enzymatic processes reliant on small, thermally-shielded pro- 
tein reaction sites and/or on local temperature gradients which 
can force cellular substrate from decoherent to coherent activity 
(cf Davies, 2004). Considering these factors, a number of sub- 
strate essential for cellular computations are already associated 
with quantum performance characteristics, such as cytoskele- 
tal lattices (Hameroff 1994; Matsuno, 2006; Craddock et al, 
2009), photosynthetic protein complexes (Hu et al., 1998; Sener 
et al., 2005), the citric acid cycle (Matsuno, 2006)and metabolism 
(Demetrius, 2003), molecular ratchets (Matsuno, 1999, 2006; 
McFadden and Al-Khalih, 1999; Patel, 2001; Cooper, 2009), 
molecule folding (Gutin et al., 1996; Cieplak and Hoang, 2003), 
synaptic boutons and vesicles (Beck and Eccles, 1992; Schwartz 
et al., 2005), long-range enzymatic activity (Frohlich, 1968, 2004; 
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see Reimers et al, 2009 for a dissenting view), odorant receptors 
(Turin, 1996; Brookes et al., 2007; Solov'yov et al., 2012), and 
second-messenger cascades (Clark, 2010a,b,c,d, 2011, 2012b). 
Quantum effects at both informational and physical degrees of 
freedom thus seem to appear in every major aspect of cell struc- 
ture and function, from sensory transduction to gene expression 
to cellular metabolism to cell motility (Clark, 2011, 2012b, c). 
However, one of many significant questions remaining to be 
answered is whether or not neurons are capable of emulating 
levels of quantum computational performance to optimize the 
fitness of cellular decisions during both normal and challenged 
cytophysiological states. Experimental and theoretical findings 
concerning aneural ciliate heuristic-guided social behaviors tan- 
talizingly imply that all eukaryotic cells equipped with cellular 
machinery for fast autocatalytic intracellular Ca^+ signaling and 
response regulation may execute quantum-efficient algorithms to 
select and implement appropriate response strategies to better 
cope with changing ambient and homeostatic conditions (Clark, 
2010a,b,c,d, 2011, 2012b, 2013a). To partly address this issue for 
neurons, I argue in the present Hypothesis and Theory article 
that intracellular store-operated Ca^^" release offers a suitable 
and common mechanism for widespread biological evolution 
and expression of Grover's quantum algorithm in cellular life. 
I begin with brief reviews of intracellular store-operated Ca^+ 
release in neurons and the basic facets of Grover's quantum algo- 
rithm. I then narrow my discussion to highlight correspondences 
between a mathematical fire-diffuse-fire model of intracellular 
store-operated Ca^+ release and Grover's quantum algorithm, 
followed by an unprecedented, if preliminary, parameterization 
of the fire -diffuse -fire model to fit Grover's quantum algorithm 
operating specifications. Lastly, I contemplate testable model pre- 
dictions and the ecological and evolutionary impact that a cellular 
version of Grover's quantum algorithm may have for healthy and 
diseased neurons and the organisms to which they are invested. 

INTRACELLULAR Ca^+ DYNAMICS AND RESPONSE 
REGULATION IN NEURONS 

Extensive varieties of functional Ca^+ channels, transporters, 
and exchangers are expressed by eukaryotic cells. Each protein 
type can be directly or indirectly involved in cellular response- 
regulatory pathways and/or Ca^+homeostasis. The main classes 
of Ca^+ channels, transporters, and exchangers across animal 
phylogeny rely on mechanosensitive [e.g., transient recep- 
tor potential (TRP) compression and stretch receptors], ATP- 
dependent [e.g., sarcoplasmic-endoplasmic-reticulum (SERCA) 
and plasma-membrane ATPase (PMCA) Ca^+ uptake/extrusion 
pumps), ion-gated (e.g., Ca^+/H+ and Na^"/Ca^+ exchangers)], 
voltage-gated (e.g., L-, N-, P/Q-, R-, and T-type receptors), 
ligand-gated [e.g., inositol 1,4,5-trisphosphate (IP3), ryanodine 
(Ry), and N-methyl-D-aspartate receptors (NMDA)], and pep- 
tidergic porin (e.g., aquaporins) mechanisms (Clark, 2013b; 
Clark and Eisenstein, 2013; Clark et al, 2013; Dolphin, 2006; 
Foskett et al., 2007). The majority of known selective or 
nonselective Ca^+ channel, transporter, and exchanger sys- 
tems are well identified and studied for a range of differen- 
tiated animal cell types, such as neurons and myocytes. With 
respect to neurons (Figure 1, left panel), scientific attention is 



frequently given to those Ca^+ -permeable or -activating pro- 
teins, such as L- and N-type voltage-gated channels, ligand-gated 
NMDA receptors (NMDARs) and a-amino-3-hydroxy-5-methyl- 
4-isoxazolepropionicacid receptors (AMPARs), and metabotropic 
glutamate receptors (GPCRs), critical for the induction and/or 
maintenance of certain forms of synaptic structural and trans- 
mission plasticity (c£ Franks and Sejnowski, 2002; Bear, 2003; 
Malenka and Bear, 2004; lino, 2006; Cohen and Greenberg, 
2008; Levitan, 2008; Yashiro and Philpot, 2008; Okamoto et al, 
2009; Catterall, 2010; Selvaraj et al, 2010; Castillo et al, 2011; 
Fioravante and Regehr, 2011; Hartmann et al., 2011; Wright 
and Vissell, 2012). However, these same and additional proteins 
are also instrumental, for example, in experience-independent 
cytokine and chemokine immunological responses, gene expres- 
sion, cellular trafficking, and homeostasis (cf Clark, 2013b; Clark 
and Eisenstein, 2013; Foskett et al, 2007), when extracellular 
Ca^^ influx and/or triggered IP3 -dependent store-operated intra- 
cellular Ca^^" -induced Ca^+ reactions (CICRs) help control endo- 
some transport, membrane remodeling, and up- and downregu- 
lation of metabolic and catabolic processes. In classic scenarios of 
facilitated and depressed glutamatergic synaptic function, respec- 
tively known as long-term potentiation (LTP) and depression 
(LTD), extracellular Ca^+ enters the post-synaptic cell through 
activated NMDARs during the induction phase of plasticity. Ca^+ 
loading in dendritic spines often evokes either LTP or LTD in a 
concentration-dependent manner, with higher and lower levels of 
Ca^+ tending to respectively produce LTP and LTD (cf Malenka 
and Bear, 2004). As free Ca^+ diffuses in spines and perhaps other 
cell compartments, various Ca^+ -dependent messenger systems 
become activated and contribute to the induction and mainte- 
nance phases of plasticity. Signal transduction by Ca^+ and its 
sensors/binding proteins, such as calmodulin and calcinurins, 
stimulate cascading enzymatic activity from calcium-calmodulin 
kinase II (CaMKII) and IV, nitric oxide synthase, protein kinase C, 
tyrosine kinase Src, mitogen-activated protein kinase, and other 
molecular complexes that lead to enduring changes through post- 
synaptic CREB-dependent transcription and immediate early 
gene activation, post-synaptic receptor synthesis, transport, and 
distribution, pre- and post-synapse geometry, and pre-synaptic 
vesicular transport and docking at transmitter release zones (cf. 
Malenka and Bear, 2004). Since the number and spatial distribu- 
tion of Ca^^" -dependent LTP and LTD events can quickly exhaust 
NMDAR transients, LTP and LTD must be supported by store- 
operated CICRs (cf Malenka and Bear, 2004; Verkhratsky, 2005). 
The initial requirements of Ca-^+ loading to excite CICRs for 
expression of LTP and LTD differ according to the frequency 
of post-synaptic stimulation, whether LTP or LTD develops, 
and the type of neuron in which they occur. However, it is 
now accepted that IP3 receptor (IP3R)-mediated CICRs assist in 
directing response regulation under physiological constraints of 
neuronal synaptic transmission and plasticity. 

IP3RS AND CYTOPHYSIOLOGY OF CICRs IN NEURONS AND 
OTHER DIFFERENTIATED ANIMAL CELLS 

Four integral 3 1 0-kDA tetrameric IP3 R isoforms, all permeable to 
certain mono- and divalent cations, especially Ca^+, are expressed 
by animals (cf Taylor et al., 2004; Foskett et al., 2007; Taylor 
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Presynaptic Cell Grover's Quantum Algorithm 




FIGURE 1 I Calcium-Induced calcium reactions (CICRs) emulate 
Grover's quantum algorithm in neuronal information processing. Left 
panel portrays major characteristic substrate (e.g., receptors, organelles, 
etc.) involved in Ca^+-mediated response regulation of arbitrary 
glutamatergic neurons, including, but not limited to, substrate critical for 
synaptic plasticity, cellular energetics, immunoprotection, homeostasis, 
gene expression, biosynthesis, molecular trafficking, cytoskeletal 
organization, and cell fate. Similar mechanisms affect both pre- and 
post-synaptic neurons, but, for descriptive purposes, post-synaptic cell 
activity is emphasized. Ca^+ entry into the post-synaptic neuron through 
voltage-gated receptor (VGC), ligand-gated receptor (LGC), and transient 
potential receptor (TRP) channels and stimulated inositol 1 ,4,5-trisphosphate 
(IP3) production by activated G-protein coupled receptors (GCR) help 
initiate cytosolic CICRs from integral IP3 receptors (IP3R) located along the 
endoplasmic reticulum (ER) membrane. CICRs may cause traveling waves 
of varying velocities and patterns which emulate search routines capable of 
eliciting/suppressing appropriate response regulation from different cellular 
compartments. Lower right panel illustrates CICR saltatory and continuous 
waves. Saltatory Ca^+ waves and the information they carry conduct at 
velocities {V) proportional to the classical Ca^+ diffusion coefficient (D). 



Whereas, faster continuous Ca^+ waves and the information they transmit 
move at velocities proportional to the square-root of the classical Ca^+ 
diffusion coefficient. Coefficient D of continuous waves for either 
intercluster or intracluster diffusion is assumed to be up to orders of 
magnitude greater than that for saltatory waves. The quadratic disparity in 
the velocities of saltatory and continuous waves corresponds to the 
root-rate increase of information processing by Grover's quantum algorithm 
over classical algorithms. Upper right panel shows schematic of Grover's 
quantum algorithm. The algorithm takes as input n qubits, upon which it 
performs Hadamard transformations (H®") and Grover's operation (GO) to 
find a target m of M solutions stored in database W. Regardless of 
whether one or more consultations of the Oracle are needed, Grover's 
quantum algorithm finds the target solution within 0= W^'^ algorithmic 
steps or operations O. Additional abbreviations: arachidonic acid (AA), Ca^+ 
binding molecule (CBM), Ca^+ uniporter (Uni), diacylgycerol (DG), Golgi 
apparatus (Golgi), L-glutamate (L-Glu), nucleus (NucI), mitochondria (Mito), 
nitric oxide (NO), nitric oxide synthase (NOS), phospholipase A2 (PLA2), 
phospholipase C (PLC), plasma-membrane Ca^+ ATPase (PMCA), ryanodine 
receptor (RyR), sarcoplasmic-endoplasmic-reticulum Ca^+ ATPase (SERCA), 
Na+/Ca^+ exchanger (Exch), synaptic vesicle (SV). 



and Tovey, 2010). Three of these isoforms, types 1 through 3, are 
encoded in vertebrate genomes (Patel et al., 1999; Taylor et al., 
1999). Invertebrate genomes encode the remaining receptor iso- 
form closely related to IP3R type 1 (Iwasaki et al., 2002; lonescu 
et al, 2006). Except for structurally and functionally similar 



polymorphic ryanodine receptors (RyRs), IP3RS are the foremost 
ion pore responsible for nonmitochondrial store-operated Ca^+ 
release in animal cells, including neurons. Receptors concentrate 
in the membrane of the ER (e.g., Ross et al., 1989; Otsu et al., 
1990; Ferreri-Jacobia et al., 2005), the principle calcium storage 
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site of most cells, but are also located in the nuclear envelope, 
Golgi apparatus, secretory vesicles, and plasma membrane of var- 
ious differentiated cell types (Ross et al, 1989; Tanimura et al, 
2000; Dellis et al., 2006). In mammalian brain, where IP3RS are 
very rich in the cerebellum (Supattapone et al., 1988; Furuichi 
et al., 1989; Sharp et al, 1999) and hippocampus (Furuichi et al., 
1993, 1994; Sharp et al., 1993), greater heterogeneity in intra- 
cellular distribution occurs for receptor isoforms. IP3R type 1, 
the most abundant isoform in brain, resides in dendrites, cell 
bodies, axons, and synaptic terminals of cerebellar Purkinje cells, 
while being largely confined to soma and proximal dendrites in 
other neurons (Ryugo et al, 1995; Dent et al, 1996; Sharp et al., 
1999). The IP3R type 3, in contrast, is localized to neuropil and 
neuronal terminals (Sharp et al., 1999). Consistent with receptor 
distributions, IP3RS act as a prominent signal interface between 
the ER and most other organelles, including mitochondria, to 
directly and indirectly affect cell processes (cf. Ponce-Dawson 
et al., 1999; Strier et al, 2003; Coombes et al, 2004; Taylor 
et al, 2004; Fraiman et al., 2006; Foskett et al, 2007; Solovey 
and Ponce-Dawson, 2010; Taylor and Tovey, 2010). Importantly, 
large observed variations in receptor structural identity may pre- 
dictably correspond to an equally large functional diversity with 
subtle developmental and physiologic consequences for specific 
IPsR-populated organs and tissues. However, receptor types likely 
provide complementary and redundant substrate for intracellular 
Ca^+ signaling. Redundancy may be apparent in small observed 
divergences in agonist binding affinities across receptor types 
(cf. Foskett et al, 2007; Clark et al, 2013). Without signifi- 
cant variance in agonist sensitivity, separate isoforms are unable 
to titrate Ca^+ -dependent physiological responses to differential 
cytosolic ligand concentrations. Moreover, channel phosphory- 
lation and protein interactions tend to cause similar respective 
aUosteric modification of heterotrophic ligand binding regardless 
of receptor type (cf Foskett et al, 2007; Clark et al., 2013). 

Well studied for oocyte development, myocardial activity, and 
cerebellar cortex neural output, IPaR-mediated ER Ca^+ release 
aids, for example, in regulating protein and lipid biosynthesis, 
cell energetics, stress responses, cell fate and death, synaptic plas- 
ticity, and immunoprotection (cf. Clark and Eisenstein, 2013). 
The complex spatiotemporal information conveyed in Ca^+ sig- 
nals is highly dependent upon IP3 production by phospholipase 
C (PLC)-P and -y, distribution and activation/inhibition proper- 
ties of IP3RS, uptake and storage of Ca^+ in the ER and other 
reservoirs, and influences over Ca^+ diffusion (cf. Clark and 
Eisenstein, 2013). The ER membrane encompasses a cisternal 
space that occupies about 10% of cell volume and importantly 
harbors cotranslational proteins, lipids, and divalent ions, such 
as Ca^+. ER-membrane-bound Ca^+ ATPases (e.g., sarcoplasmic- 
endoplasmic-reticulum Ca^+ ATPase) sequester free cytosolic 
Ca^+ in the ER lumen, where it stays free or becomes attached 
to buffers. Estimates of total luminal Ca^+ concentrations are 
as high as 1 mM. The fraction of unbound luminal Ca^^ ranges 
from 100 to 700 (iM (Montero et al, 1995; Bygrave and Benedetti, 
1996; Pinton et al, 1998; Alvarez and Montero, 2002; Bassik 
et al, 2004; Palmer et al, 2004; Verkhratsky, 2005). Although 
ER Ca^"*" modulates IP3R function from the luminal side via 
interactions with channel sensor domains, cobinding of IP3 



and Ca^+ to cytosolic sites opens the gated IP3R channel in 
a concentration-dependent manner, driving passive Ca^^ flux 
down its electrochemical gradient and into the cytosol. In IP3 
concentrations ranging between 100 nM and 1 (xM that contin- 
uously saturate receptor binding sites for IP3, Ca^+ generally 
excites and blocks IP3R activity in respective low (e.g., ~50nM- 
1 |xM) and high (e.g., >10 |xM) concentrations (cf. Foskett et al., 
2007). Changes in IP3R conformation and pore permeability 
occur due to IP3 and Ca^+ aUosteric interactions that disso- 
ciate suppressor, calmodulin, and gatekeeper receptor regions, 
repositioning the transmembrane gate and activating Ca^+ con- 
ductance (cf. Foskett et al., 2007; Clark and Eisenstein, 2013). In 
absence of IP3 binding, low-affinity binding of cytosolic Ca^+ 
to one of two calmodulin heads occludes the ion channel as 
calmodulin crosslinks with suppressor and gatekeeper regions 
of adjacent receptor subunits. The receptor lumen stays closed 
and inactive when only cytosolic IP3 binds to receptor sites. 
Depending on recording preparations, receptor type, and other 
factors, Ca^+ conductance and current through single channels 
have been electrophysiologically measured at around 10-125 pS 
and 0. 1-0.5 pA, respectively, (cf. Foskett et al., 2007). Maximum 
mean duration of IP3R opening tends to be no more than 15 and 
40 ms for respective vertebrate and invertebrate receptor types (cf. 
Foskett et al, 2007). The lengths of these periods are indepen- 
dent of agonist concentrations. But the durations of subsequent 
prolonged refractory or reversible inactivation periods are deter- 
mined by agonist concentrations. Frequency of channel activity 
and graded Ca^+ mobilization are thus primarily due to cyto- 
plasmic ligand levels, with additional regulation by nucleotides, 
phosphorylation, redox states, and protein interactions. 

Stimulation of an individual or several IP3RS evokes a spa- 
tially discrete Ca^+ release usually termed a blip (Figure 1, lower 
right panel). Blips are the most elemental Ca^^" release event (cf. 
Foskett et al, 2007; Solovey and Ponce-Dawson, 2010). Quantal 
release creates a microdomain of high cytosolic Ca^+ concentra- 
tion that can exceed 100 |a,M near the opening of an activated 
channel (Naraghi and Neher, 1997; Rios and Stern, 1997; Neher, 
1998). Once a 10-|xM Ca^+ threshold is reached, free cytosolic 
Ca^+ begins to exert inhibitory feedback control over all-or-none 
openings of surrounding IP3RS. Inhibitory control of IP3RS is 
proportionally tuned by presence of local IP3 concentrations, 
which interfere with the ability of Ca^+ to bind to low-affinity 
sites. Furthermore, rapid buffering by both mobile and immobile 
Ca^+ traps often limit diffusion of free cytosolic Ca^^ to a radius 
no greater than 5 [xm from the source channel (Allbritton et al., 
1992). Buffers alone are typically insufficient to quench a Ca^+ 
release event involving multiple adjacent receptors. But beyond 
the distance of 5 yim, cytosolic Ca^+ concentrations, homeostat- 
ically regulated by transporters, exchangers, and porins, steeply 
decline from 1 |xM to ~50 nM without widespread Ca^^^ mobi- 
lization (Naraghi and Neher, 1997; Rios and Stern, 1997; Neher, 
1998). Since IP3RS organize in autocatalyst-linked clusters on the 
ER membrane, a larger coordinated Ca^+ release event, com- 
monly called a puff or spark, can be initiated following a blip 
(cf. Ponce-Dawson et al, 1999; Strier et al, 2003; Coombes et al., 
2004; Fraiman et al., 2006; Foskett et al, 2007; Solovey and 
Ponce-Dawson, 2010). Puffs occur as Ca^+ from the blip diffuses 
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and binds to neighboring inactive IP3RS already docked by IP3 
(Figure 1, lower right panel), an autocatalytic event referred 
to as a CICR or fire-diffuse-fire reaction (cf Yao et al., 1995; 
Ponce-Dawson et al, 1999; Strier et al., 2003; Coombes et al., 
2004; Taylor et al, 2004; Fraiman et al., 2006; Guisoni and de 
Oliveira, 2006; Shuai et al., 2006; Foskett et al, 2007; Bruno 
et al., 2009; Smith and Parker, 2009; Ur-Rahman et al., 2009; 
Solovey and Ponce-Dawson, 2010; Taylor andTovey, 2010). Single 
IP3RS are typically separated by 10-20 nm within a cluster of 
perhaps 50 receptors covering a maximum ER-membrane sur- 
face area of 400 nm^ in some cells (Shuai et al, 2006; Bruno 
et al., 2009; Ur-Rahman et al., 2009). Individual clusters can 
be separated by regular or irregular distances of up to around 
2 (Jim (Yao et al., 1995). Though cytosolic Ca^+ binding proteins 
may alter Ca^^" transport in the vicinity of a receptor chan- 
nel cluster, the interchannel distances are too short to prevent 
most ions from diffusing. Therefore, when an estimated 20 to 
35 IP3RS become simultaneously bound with coligands, Ca^+ 
puffs arise (Shuai et al., 2006; Bruno et al., 2009; Smith and 
Parker, 2009). As mentioned earlier, this autocatalytic coupling 
of clustered receptors forms local and global cellular networks 
or lattices that can generate either small- or large-scale mobi- 
lization of Ca^+. Irregularities in spatial organization of IP3R 
clusters together with channel coupling associated with CICRs 
produce assorted intracellular Ca^+ signals. The amplitude, fre- 
quency, and velocity of signals vary according to cytosolic Ca^+ 
buffer concentrations, feedforward excitation of IP3RS via possi- 
ble high affinity Ca^+ binding, feedback inhibition of IP3RS via 
possible low affinity Ca^^ binding, and crosstalk with additional 
messenger systems (cf Clark, 2011, 2012b), such as cAMP path- 
ways (e.g., Siso-Nadal et al, 2009). Puffs may trigger global Ca^+ 
waves as cellular conditions promote sustained CICRs. Waves can 
oscillate, extinguish, and travel throughout different cell com- 
partments as saltatory, continuous, or anisotropic fronts. The 
spatiotemporal variability of blips, puffs, and large-scale waves 
suggests a high degree of specificity is achieved for intracellular 
Ca^^ signaling, reducing the likelihood of corruption and loss 
of transmitted information content by noisy intracellular pro- 
cesses (cf Clark, 2011, 2012b, 2013a). In many respects then, 
the behavior of neuronal CICRs conforms to fundamental prin- 
ciples and attributes of (classical and quantum) search algorithms 
and patterns used to efficiently find and execute various kinds of 
appropriate cellular responses to extracellular and/or intracellu- 
lar stimuli (Clark, 2010a,b,c,d, 2011, 2012a,b, 2013a; Clark and 
Eisenstein, 2013). 

DESCRIPTION OF GROVER'S QUANTUM ALGORITHM 

Before entering into discussion on the relationship between 
CICRs and Grover's quantum algorithm, I now identify basic 
specifications of Grover's quantum algorithm through a short 
primer of information and computational theory. In standard 
quantum information theory and computation, the classical bit 
originated by Shannon (1948a,b) is replaced with the quantum 
bit or qubit, a concept and term, respectively credited to Weisner 
(1983) and Schumacher (1995). Qubits are information units that 
maybe transmitted, transformed, stored, and measured. The pos- 
sible states of a single qubit in Dirac notation are the orthonornal 



"1" 




"0" 


and |1) = 


_0_ 


_1_ 



unit vectors or basis (eigen)states |0) 

which span a two dimensional vector, state, or Hilbert space. 
States |0> and |1> correspond to classical bit states of 0 and 1. 
However, unlike classical bits, these states may form an indefinite 
linear combination or superposition: = a\0) + b\l), where 
variables a and b are complex numbers called vector ampli- 
tudes (cf Nielsen and Chuang, 2000). By way of a common 
example useful to later exposition of Grover's quantum algo- 
rithm, states |0) and 1 1) forming the unique quantum supposition 
(|0) - |l))/2'''^ have respective amplitudes 1/2^/^ and-l/2^/^If 
qubit states remain linearly independent, then any measurement 
or eigenvalue m on the system | i/f,> via arbitrary Hermitian mea- 
surement operators M," = andMo = /— Xii^to \i^t)(i^i\> 

1 0" 
0 1 



where I is the identity matrix 



will decompose into a sin- 



gle pure state |0) or |1> with respective probabilities \a\^ and 
\b\^, so that measurement probability p,- = (Vf,|M, = \a\^ + 
\b\^ = 1 (cf Nielsen and Chuang, 2000). Conversely, in agree- 
ment with the entropic uncertainty principle, an informational 
analog to Heisenberg's uncertainty principle, indistinguishable 
or nonorthonormal quantum states cannot be measured with 
certitude because of variances inherent in observables. As one 
may expect, these concepts become slightly more complicated in 
instances of composite or joint systems — those involving two or 
more qubits. Although the full implications of composite systems 
go beyond the scope of this article, it is important to intro- 
duce some content on the subject for future consideration. Take 
two entangled component systems, each respectively described 
by superposition states |i/f,) = (|0) -|- |l))/2'''^ and \ = (|0) - 
|l))/2^/^. The state space of such a bipartite system is defined 
by the tensor product \ (gi = \ = (|00) - |ll))/2'/^ 
yielding, in the present case, the third Bell basis state or Einstein- 
Podolsky-Rosen pair. Composite quantum information systems, 
such as the four Bell basis states or Einstein-Podolsky-Rosen 
pairs, play significant roles in superdense coding, information 
encryption, error diagnosis and correction, and other aspects of 
quantum computation, including execution of quantum algo- 
rithms (cf Nielsen and Chuang, 2000). 

As previously noted, quantum computational methods entail 
use of specialized quantum gates and circuits to form algo- 
rithms that manipulate qubits to purposefully arrive at some goal 
state, much as would be accomplished for classical bits operated 
on by classical logic gates and circuits. Purposes may include, 
for instance, finding correct solutions to difficult or classically 
intractable factorization, ordering, counting, and search prob- 
lems. Now imagine a quintessential large database or map of 
salient landmarks, such as cities (or, as will be detailed below, 
rate-limiting parameters for selective spatiotemporal chemical 
diffusion patterns). Using a classical search algorithm to dis- 
cover the shortest route n among all possible routes N through 
every city on the map, a dilemma known as the Hamiltonian 
cycle decision problem, requires N total operations O or algo- 
rithmic steps [i.e., 0(N)]. The same search problem may be 
accelerated to 0{N^^^) with Grover's quantum algorithm and 
its unique Grover's operator (Grover, 1996) (Figure 1, upper 
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right panel). Typical uses of Grover's quantum algorithm search 
register entries i = {;„!« = 0, . . . , N — 1} indexed to actual ele- 
ments n of N. This convention is created so that the database 
can be conveniently set to N = 2" bits of storage, with a sub- 
set M = 1 < M < iV = {im|m = 1 < m < N — 1} of exact solu- 
tions. Another convention enlists a function 9 that accepts as 
input an integer x = i valued over the range 0 to N — 1 . The 
function returns results 9{x) = 1 or 6{x) = 0 when a solution 
m is or is not obtained, respectively. Given these constraints, 
Grover's operator consists of four distinct sequential procedures 
(cf. Nielsen and Chuang, 2000): (1) application of the Oracle, (2) 
application of the Hadamard transformation, (3) application of 
a conditional phase shift, and (4) application of the Hadamard 
transformation. The operator first samples inputs from the ini- 
tial equally weighted superposition state = 11x=o I*) 
generated by the efficient Hadamard transform, H®", then it 
labels problem solutions through the unitary action of the 
Oracle: 

\x)im - |l»/2l/2 ^ -l'*"'W(|0> - |l»/2^/^ (1) 

where \x) represents the index-register qubit set to 
|0) and (|0) - |l>)/2'/2 is the single Oracle qubit which 
assists in flipping or phase-shifting \x) only when 9ix) returns 
1 as a result. A second Hadamard transformation is utilized fol- 
lowing the Oracle call to place labeled qubits into superposition. 
From this state, a conditional phase shift, \x) —(—1)^'^^ \x), 
becomes executed for all basis states not equal to |0). The final 
Hadamard transformation again puts the register qubit into 
superposition for possible further Oracle summons in the event a 
target solution is not located, although successful search attempts 
may require no more than one Oracle call. 

CORRESPONDENCE BETWEEN CICRs AND GROVER'S 
QUANTUM ALGORITHM 

If modifiable operation of linked intracellular Ca^+ release sites 
and associated affector/effector systems function as quantum-like 
computational networks for response regulation, as reported by 
Clark (201 Oa,b,c,d, 20 11 , 201 2a,b, 20 1 3a) , then reaction-diffusion 
equations should reveal classical and quantum properties of 
search algorithms and of search patterns selectively applied to 
those networks. One such reaction-diffusion equation, the simple 
fire-diffuse-fire model of Ca^+ propagation, defines Ca^+ waves 
by the following evolution equation (Ponce-Dawson et al., 1999): 

3 [Ca^+] /dt = D{d^ [Ca^+j /d^x) + {a /(fx) 

^S(x- Xi)H{t - ti)H{t, + r - t) , (2) 

where [Ca^^](x, f) is the average concentration of calcium in 
directions perpendicular to the direction x of propagation, S{^) 
is the S function, H(f ) is the Heavyside step function (and not 
the Hadamard transform), D is the classical Ca^+ diffusion coef- 
ficient, d is the (mean) distance between Ca^+ release sites, t, is 
the first time the ith Ca^+ release site reaches threshold value, a 
is the total number of Ca^+ ions released per storage site, and r is 
the period that receptors remain open to release Ca^+ in a single 



event. Notably, the evolution equation and parameters for this 
sort of model are general enough to fit conditions stipulating IP3R 
clusters or individual IP3RS as Ca^^ release sites (cf. Guisoni and 
de Oliveira, 2006; Solovey and Ponce-Dawson, 2010), with corre- 
sponding changes in intersite distances, timescales, and released 
Ca^+ concentrations. The dynamics of the equation additionally 
rely on two dimensionless parameters, F and /J (Ponce-Dawson 
et al, 1999), with: 

r = (a/d^)/([Ca^+]^-[Ca^+],), (3) 

where a/cP is the concentration of Ca^+ released and [Ca^+]r 
and [Ca^+];, are respective threshold and basal concentrations of 
Ca^+, and 

l3 = {D/r)/d\ (4) 

Parameter P of Equation 3 acts as a multiplicative variable gov- 
erning the ease of starting a Ca^+ wave, the velocity at which 
it will travel, and, therefore, the effectiveness of the wave to 
effect response regulation. Whereas, the value of parameter f} of 
Equation 4 defines whether a Ca^+ wave propagates with a slow 
saltatory, fast continuous, or intermediate mixed front. As chemi- 
cal processes, such as the duration of Ca^+ mobilization or degree 
of Ca^+ buffer overload, become rate-limiting over interstore 
distances, 5S> 1 and Ca^+ waves transition from saltatory to 
continuous waves (Figure 1, lower right panel). Moreover, salta- 
tory wave propagation travels at a rate proportional to the Ca^+ 
diffusion coefficient, v~ {D/d)g^^r, where g^^ is an inverse 
function (Ponce-Dawson et al, 1999). In contrast, the velocity of 
continuous waves is proportional to the square-root of the Ca^+ 
diffusion coefficient, v ~ (D/r)^^^f^^r, where /^^ is an inverse 
function (Ponce-Dawson et al, 1999). When continuous waves 
have sufficiently large F, wave velocity approximates the Luther 
equation, v = a(D/r)^^^ with a = F^^^. The equations for wave 
velocity should not be misleading, as continuous waves are often 
faster than saltatory ones (cf Izu et al., 2001). (Noted excep- 
tions include the large-scale completely homogenous fertilization 
waves of oocytes.) If the value of D is the same for both saltatory 
and continuous waves and t is much larger for continuous waves, 
then saltatory Ca^+ waves would always transmit at faster speeds. 
But when compared to saltatory waves, continuous waves well 
exceed buffering capacities of slow and fast Ca^^ -buffer species 
and display far greater diffusion coefficients and shorter intersite 
diffusion times (cf Strier et al., 2003) for both intercluster and 
intracluster models of diffusion for constant small d. The effect 
of overcoming buffering capacity on continuous wave velocity 
may be also amplified by diminution of r on the order of one 
to two magnitudes to the ms timescale (cf Izu et al., 2001; Strier 
et al, 2003; Foskett et al, 2007), which differs from r given by 
Ponce-Dawson et al. (1999). In either situation of saltatory or 
continuous waves, wave conduction generally fails for small D 
(e.g., <10 [irn^/s), large d (e.g., d > 3 |xm), and extremely small 
or large F (cf. Keizer et al., 1998; Ponce-Dawson et al, 1999; 
Strier et al, 2003). These values reflect significant differences in 
the physiological roles of saltatory and continuous Ca^+ waves 
(e.g., Keizer et al., 1998), with the former believed to inhibit local 
and global cellular responses via Ca^+-wave conduction failure 
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and CICR blockade and the latter believed to evoke and integrate 
a range of local and global cellular responses through complex 
spatiotemporal patterns and widespread delivery of information 
throughout the cell. Fire-diffuse-fire models of Ca^+ propagation 
are remarkably robust, explaining the diffusive characteristics of 
store-operated Ca^+ regulation in a generous variety of eukary- 
otic cell types, including oocytes, cardiac myocytes, and neurons 
{e.g., Ponce-Dawson et al, 1999; Strier et al, 2003; Coombes 
et al, 2004; Timofeva and Coombes, 2004; Fraiman et al., 2006; 
Guisoni and de Oliveira, 2006; Thul et al., 2007; Solovey and 
Ponce-Dawson, 2010; Bressloff 2014). 

Importantly, the leading edge of intracellular Ca^+ gradients 
passing between separate receptor clusters is relatively slow, being 
experimentally recorded to typically travel at velocities rang- 
ing from 20 to 70|im/s (laffe, 1993). Some reports, however, 
indicate continuous longitudinal waves can approach an aston- 
ishing 6000 |im/s in live cells (Miura et al., 1999). In any event, 
when N and D are numerically related (Clark, 2010a, 2012b), 
the quadratic disparity between the speeds (i.e., intersite-travel 
time) of saltatory and continuous Ca^+ waves appears consis- 
tent with expression of a square-root quantum algorithm that 
increases target searches, such as searching for the most appro- 
priate response to external and/or internal stimuli, in living cells 
beyond the bounds of classical algorithms (Clark, 2010a,b, 2011, 
2012b) (Figure 1, lower and upper right panels). In the algorith- 
mic (rather than strictly physicochemical) sense, the properties 
of classical Ca^+ diffusion effectively describe a phenomenolog- 
ical basis for Grover's quantum algorithm (cf Clark, 2010a,b). 
A physicochemical manifestation of quantum mechanics via a 
quantum diffusion term is unnecessary to produce quantum- 
efficient algorithm searches, a result perhaps counterintuitive 
for most biophysicists since diffusive processes may be classical, 
quantum, or semiclassical/semiquantum in nature. At biologi- 
cally relevant subsecond times, warm temperatures, and microm- 
eter scales described for the conditions of Ca^+ fire-diffuse-fire 
reactions (cf, Ponce-Dawson et al., 1999), quantum diffusive 
processes produce minor effects without thermodynamic shield- 
ing, such as in the case of bacterial photosynthetic reaction 
cores (Hu et al., 1998; Sener et al, 2005), or pump-process 
energy transfer, such as in the possible case of actomyosin poly- 
merization (Matsuno, 1999). Neither thermodynamic shielding 
nor pump-process constraints must occur for initiation and 
maintenance of classical fire-diffuse-fire reactions. The diffusion 
coefficient or diffusivity, D, in Equations 2 and 4 is a purely 
classical parameter generally derived from Pick's laws and the 
Einstein-Smoluchowski relation as D = /x/cbT, where fi is particle 
mobility or the inverse drag coefficient, is Boltzmann's con- 
stant, and T is temperature in degrees Kelvin (cf. Clark, 2012b). 
This equation, via the Einstein-Sutherland relation, becomes the 
Einstein-Stokes equation in one dimension, D = k^T /In-qr, and 
in three dimensions, D = k^T/6ji7ir, for spherical particles of 
radius r moving through a fluid of viscosity at a low Reynolds 
number. If the density of the diffusing material affects D, then 
the diffusion equation is nonlinear and D is taken to be vari- 
able (cf Clark, 2012b). When independent of thermodynamic 
influences, D = h/Ajtm, where h is Planck's constant and m is 
the mass of the diffusing particle. The latter definition of the 



diffusion coefficient is quantum mechanical (cf Clark, 2012b). 
One may draw direct comparison of these sorts of effects with 
the operation of closely related technological quantum networks 
performing search functions (e.g., Bianconi and Barabasi, 2001; 
Bianconi, 2002a,b, 2003; Stella et al, 2005; Clark, 2010b,c,d, 
2011, 2012a,b, 2013a). In such instances, observed statistical 
quantum-like outcomes, often referred to as quantum mechanical 
analogs, emerge from weighted macroscale computational net- 
works and their parameters capable of both classical and quantum 
behavior. For technological systems (Bennett, 2003; Ladyman 
et al, 2007) and individual cells (Clark, 2010b,c,d, 2011, 2012b, 
2013a; Berut et al, 2012; Mehta and Schwab, 2012), this behav- 
ior is consistent with Landauer's principle of energy/information 
transfer. Similar to computational network analogs of quan- 
tum behavior, the classical Ca^+ diffusion coefficient might 
instantiate a computational analog of quantum mechanical sys- 
tems without actually residing in a physicochemical quantum 
regime. 

Although the physical expression of the diffusion coefficient 
should be entertained, it must he stressed that in some sense any 
reflection is superfluous with respect to application of Grover's 
quantum algorithm. The reason for this, as indicated in the pre- 
vious paragraph, is that a relationship between N and D seems 
apparent (and wOl be established in below sections). The vari- 
able N represents the total number of search elements queried 
by Grover's quantum algorithm. Its value is neither quantum nor 
classical! Therefore, D needs to be neither quantum nor classi- 
cal to effect Grover's quantum algorithm in a cellular system! 
That is, what makes Grover's quantum algorithm quantum in 
nature is its action on a search field, not necessarily the prop- 
erties of the search field itself In view that classical diffusion 
terms satisfy the quadratic improvements needed for Grover's 
quantum algorithm, it is interesting that search selectivity by 
such an algorithm in single cells may be enhanced, instead of 
being damped, by diffusion barriers sometimes causing unsta- 
ble Ca^+ gradients. Anisotropic patterns of Ca^^" diffusion due 
to free cytosolic chaperons and buffers (Chen et al., 2008, 2009) 
and frequency and amplitude modulated Ca^^" liberation (De 
Pitta et al., 2008, 2009) have been reported to help improve the 
specificity of encoding sensory information transmitted by intra- 
cellular Ca^+ cascades. Intracellular spaces are filled with Ca^+ 
traps, such as immobile binding sites. Though traps reduce the 
effective diffusion coefficient below expected values for free dif- 
fusion, propagation of information can move faster than single 
particle diffusion (Pando et al., 2006). The storage and retrieval 
of that information is expected to be further refined by recur- 
sive phosphorelays affecting Ca^+ permeability to extracellular 
sources and subsequent reactivation CICRs (Clark, 2010a,b,c,d, 
2011, 2012a,b, 2013a). Collectively, these findings should be put 
into local and global contexts of mixed and continuous wave 
fronts rather than taken to mean that spatiotemporal patterns 
of ineffectual punctuate store-operated Ca^+ emissions or unre- 
liable saltatory Ca^+ waves serve as useful media to convey 
information vital to response regulation (cf Keizer et al., 1998). 
These events likely rather function as wave guides that direct 
mixed or continuous waves to specific target locations within cell 
compartments. 
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PARAMETERIZING THE NEURONAL FIRE-DIFFUSE-FIRE 
MODEL FOR GROVER'S QUANTUM ALGORITHM 

The superficial similarity between variables N and D in quadratic 
processing efficiency encourages mathematical treatment that 
supplants analogy and parameterizes the Ca^+ fire-difi^use-fire 
model to more precisely fit specifications of Grover's quantum 
algorithm (Clark, 2011, 2012b). Because the model captures local 
dynamics of individual IP3R permeability, which inherently con- 
trol the evolution of collective wave behavior within and across 
networked IP3R clusters, Grover's quantum algorithm may be 
reduced to the scale of a single receptor channel, as opposed 
to receptor clusters or an entire ER membrane of receptors. 
This convention has several attractive qualities. First, an inabil- 
ity to make quantum measurements due to quick decoherence 
rates of superposed states becomes a negligible confound for 
quantum logic operators the physical size of tetrameric protein 
channels (cf Beck and Eccles, 1992; Gutin et al., 1996; Turin, 
1996; Cieplak and Hoang, 2003; Davies, 2004; Brookes et al, 
2007; Solov'yov et al., 2012). Second, uncovering (or framing) a 
reasonable search goal and associated parameters becomes more 
conceptually manageable. In regard to both issues, the natural 
inclination would be to perhaps equate N from the previous 
example of the Hamiltonian cycle decision problem to the total 
number of possible spatiotemporal patterns of chemical diffu- 
sion needed to evoke a proper IPaR-mediated neuronal response 
to external and/or internal perturbation, with the target solution 
being the shortest chain or route of networked receptor clusters 
across the entire or a circumscribed area of the ER membrane 
surface. While this decision-problem situation correctly assumes 
that some or all receptor clusters can be activated as a collec- 
tive search routine to evoke an arbitrary desired cell response, 
it challenges the spatiotemporal limits of superposed (physi- 
cal and not network-analog) quantum states needed to execute 
Grover's quantum algorithm, since each state embodies a long- 
range networked pattern of catalyst-linked receptor clusters (e.g., 
minimally greater than 2 or 4 [im) and, consequently, probably 
surpasses Wigner's mass-time uncertainty for periodic mechani- 
cal processes (cf Wigner, 1957, 1981; Reimers et al, 2009). The 
mechanism and processing efficiency of this kind of emergent 
search algorithm likely would be characteristically (physicochem- 
ically) classical in nature. Whereas, the almost instantaneous and 
continuous chemical diffusion in an approximate 10- to 20- 
nm distance between two activation-primed IP3RS (Shuai et al., 
2006; Bruno et al, 2009; Smith and Parker, 2009) avoids violat- 
ing Wigner-type quantum boundaries (cf. Pesic, 1993; Schwartz 
et al., 2005) and yields suitable conditions for expression of a 
quantum algorithm. Moreover, in terms of second messenger 
reaction-diffusion cascades and response regulation, stochastic 
punctate intracluster spatial patterns of diffusion often play sig- 
nificant roles in spark and wave initiation and are complemented 
by intermediate and overloaded ion-concentration magnitudes, 
which help force activation gradients, fast conduction velocities, 
and stable propagation to initiate cell responses (cf Keizer et al., 
1998; Izu et al., 2001; Strier et al, 2003; Chen et al, 2008, 2009; 
Solovey and Ponce-Dawson, 2010). Stipulating individual IP3RS, 
and their quantum-mechanical small reaction sites and gating 
kinetics (cf Ahern et al, 2009; Chan et al, 2010; Pitt et al, 2010; 



Li et al., 2013), as the computational apparatus of Grover's quan- 
tum algorithm allows one then to resolve the search problem to 
finding the shortest time or fastest rate m taken to reach and 
autocatalytically activate a neighboring receptor from among a 
database N of all possible time and activation outcomes (Table 1). 
Given slower intracluster saltatory Ca^+ waves have high prop- 
agation and response-regulation failure rates (e.g., Guisoni and 
de Oliveira, 2006; Solovey and Ponce-Dawson, 2010), target solu- 
tions wOl be elements of the set M of fast intracluster continuous 
Ca2+ waves. Such a query is idealized by the search for the max- 
imum Ca^+ diffusion coefficient Dmaxi which, in integer form 
spanning the range of possible integer diffusion coefficients, bears 
likeness to N. However, as individual IP3RS do not actually detect 
diffusion coefficients, the search must be conducted over a "reg- 
ister" indexing each different D with a concentration-dependent 
parameter biologically associated with D and germane to IP3R 
reaction kinetics. 

To make the present exposition more explicit and amenable 
with the previous coverage of Grover's quantum algorithm, 
I now mathematically define critical variables N, M, and D 
(Table 1), giving fuller attention to physical descriptions of the 
Grover's quantum algorithm workspaces and operators in sub- 
sequent paragraphs (Figure 1, upper right panel). As before, 
N = = 0, . . . , 00}, but, as will be determined below, is 
realistically a finite interval with respect to reaction-diffusion 
parameters. Elements «j of N are indexed with register entries 
= [i„\n = 0, . . . , N - 1}, with a subset M=1<M<N = 
[i^\m = 1 < m < N — 1} of exact solutions to the search prob- 
lem. Algebraic manipulation of Equation 4 yields the general 
definition of the classical Ca^+ diffusion coefficient: 



D = (pd^) /r. 



(5) 



However, fluctuations in loading of cytosolic Ca^+ buffering 
systems and, more appreciably, choice of rate-limiting param- 
eters (inter-IP3R) cf, r, and a for the fire-diffuse-fire model 
make the diffusion coefficient variable. Sharper definitions of D 
may be obtained from the wave velocity equations for saltatory 
Ca^+ waves, v = (D/d)g^^r, and continuous Ca^+ waves, v = 
(D/T)^/^/^^r. It follows that the respective diffusion coefficients 
for saltatory and continuous Ca^+ 



waves are: 



Ds = vdg/ r = {Dsi\l = min, . . . , max}, 
with Ds^ > 0 and Ds„„ « Dc^^, and 

Dc = r {vf/rf = {Dq\1 = min, . . . , max} 



(6) 



(7) 



with Dc„i„ » 0, Ds^,, and Dc^^ « 00. Setting = 0 < 

Dm3.^<Dc^^ = {Dnj\j=0,...,Dc^} and considering the 
one-to-one discrete mapping D^ax onto N, Dm^x —^N,M 
therefore transforms into the subset Dc = Ds^^ <^ M < N = 
{im\m = Ds^^^ <^ m < N — 1} of all correct solutions involving 
only intracluster continuous Ca^+ wave fronts. A major result 
from this interpretation is that, in order to arrive at a solu- 
tion m = Dc, quadratic and exponential speed-ups in respective 
algorithmic search time and wave velocity must coexist, with a 
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Table 1 | Comparison of major Ca^'*' fire-diffuse-fire model and Grover's quantum algorithm parameters. 

Ca^"*" fire-diffuse-fIre model* Grover's quantum algorithm 

DIFFUSION ^^^^^HH^^^^^^^^^^^^^^^^^^^^^^^^^^^II^^I^^^^ 

D={l}d^)/z = fikeT 

Ds = vdg/r = [Ds,\l= min, .... max), with Ds^,„ > 0, Ds^„ « Dc^,^ 
Dc = rivf/ rf = {Dq\I = min. . . . , max), with Dc^.„ » 0, Ds„,, and Dc„, « oo 

Drnax = 0 < Dmax < Dc„,, = {Dr,^\j = 0, . . . , Dc^^J 



WAVE DYNAMICS^^^^H 

P = (D/r)/d2 

= (Ds/T)/d2 
/8c = (Dc/T)/d2 
ySmax = (Dmax/rj/d^ 

r = (a/c/3)/([Ca2+]7--[Ca2+]b) 
rs = Tvg/^d={rs,|/=min, ... 
Tc = vf/p^'^d = jrc,]/ = min, . , 

Tmax = j/n|n= 0, .. . , Tc™,) 



max), with Tg^ 
. , max), with rr:„ » 0, rs^„ and Tc^,, « oo 



Dc = (Ds,3, «M<N) = = Dg,,, « m = W - 1 } 
Dmax ^ W = {nyiy =0, . . . , Dc^^J, one-to-one discrete 
nearest integer mapping 



rs = /F = |/f|f = 0<f<rc^J 

rmax = /=|/n|n=0. ...,W-1) 



'Model variables are suitable for describing both local (interreceptor or Intracluster) and global (Intercluster) Intracellular waves. 



maximum algorithmic search time of 0(Dmax/-Dc)^^^ and wave 



velocity of vc : 



In addition, Equations 6 and 7 are espe- 



cially useful for purposes of implementing Grover's quantum 
algorithm, since they are inversely proportional to dimensionless 
parameter F rather than p. As shown in Equation 3, reliance upon 
r permits Ca^+ diffusion coefficients to be indexed to physiolog- 
ically pertinent ratios involving Ca^+ concentrations liberated by 
IPaR-dependent store operation and free cytosolic Ca^+ concen- 
trations sensed by multiaffinity IP3R cytosolic Ca^+ binding sites. 
So that, index values of register entries i corresponding to Dm^x 
and Dc may be practically redefined by substituting Equation 5 
into Equations 6 and 7, then solving for T for separate saltatory 
and continuous Ca^+ wave forms: 



Ts = T:vg/Pd = {rs,\l = min, . . . , max} 



(8) 



with Vs^,, > 0 and r^,, « Vc^,,, and 

Tc = vflP^l^d = {FqII = min, . . . , max}, (9) 
with Tr. » 0, rs^,„ and < 00. Therefore, index 



i = {in\n = 0, 



,N-l} = rn,ax = {i„|n = 0, 



and 



solution subset index ijvf = {im\m = Tc^.^ <m< Tc^^}. 
Remember also from earlier reviewed content that Grover's 
quantum algorithm employs a function 0 that accepts as 
input an integer x = i valued over the range 0 to N — 1. The 
function returns results 6(x) = 1 or 9(x) = 0 when a solution 
m is or is not acquired, respectively. In parameterizing the 
fire-diffuse-fire model for Grover's quantum algorithm, elements 
of the index set i yielding 0{x) = 1 readily signify the solution 
subset iM = {im\m = rc^.^<m<rCr^^} for intracluster 
continuous Ca^+ waves, whereas elements of the index set i 
yielding 6{x) = 0 signify the incorrect-solution or failure subset 
ip = {ijr\f = 0 < f < Fc^j^} for aU intracluster noncontinuous 
Ca^+ waves, including saltatory and possibly mixed Ca^+ wave 
fronts. 



With key fire-diffuse-fire model parameters written in terms 
of Grover's quantum algorithm, candidate physicochemical sub- 
strate for algorithm operators can be next identified to a first 
approximation using known IP3R molecular biology and func- 
tion (Figure 2). Grover's quantum algorithm, as previously men- 
tioned, requires five distinct sequential procedures. The first of 
these steps, in notation consistent with the fire-diffuse-fire model, 
is application of the Hadamard transformation, H®"', which pre- 
pares the algorithm in an initial equally weighted superposition 

state |iA> = l/J^max'^^I]f.^o "' I*)- S^'^^ ^ ^^^^^ presumes that 
a single IP3R acting as Grover's quantum algorithm is capable 
of simultaneously detecting any probable index values Fj^nx = 
li„} marking Dmax = {D„j} and associated with x, the variable 
denoting a successful or unsuccessful search for shortest times 
or fastest rates m needed for store-released concentrations of 
free cytosolic Ca^+ to continuously diffuse and autocatalyti- 
caUy activate a nearest neighbor receptor and, thereby, ensure 
fast cellular response regulation. Furthermore, IP3R molecular 
regions and all other substrate effecting the superposition interro- 
gate the superposition state with interaction-free measurements 
to maintain the superposition state until a solution is deter- 
mined (e.g., Hosten et al., 2006). This computational feat is, of 
course, accomplished via inferential measurement of Dmax via 
measurement of index Fmax inherent in the operation of IP3RS. 
Although IP3RS may switch between four different conforma- 
tional states, the inactivated IPs-bound IP3R conformation is 
perfectly suited for the initial superposition state of Grover's 
quantum algorithm because all saltatory and continuous waves 
affiliated with Fmax retain equal probabilities of manifesting 
(Figure!). 

After initializing the IP3R into this IPs-saturated superposition 
state, which may remain indefinitely so in saturating IP3 con- 
centrations, the final four steps of Grover's quantum algorithm 
involve execution of Grover's operator — application of the Oracle, 
a second Hadamard transformation, a conditional phase shift, 
and the last Hadamard transformation. Grover's operator may 
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FIGURE 2 I Model of conformation, ion permeability, and corresponding 
Grover's quantum-algorithm function of an inositol 1,4,5-trisphosphate 
receptor channel (IP3R). Each cross-section contains two of four complete 
IP3R subunits. When only cytosolic IP3 (blue sphere) binds, the receptor 
lumen stays closed and inactive. Cobinding of cytosolic IP3 and Ca^+ (red 
sphere) to separate high-affinity sites proximal to the IPs-binding domain 
dissociates suppressor (Suppr), calmodulin (CaM), and gatekeeper (Keeper) 
regions, repositioning the transmembrane gate (Gate) and activating Ca^+ 
conductance. In absence of IP3 binding, low-affinity binding of cytosolic Ca^+ 
to one of two calmodulin heads occludes the ion channel as calmodulin 
crosslinks with suppressor and gatekeeper regions of adjacent receptor 
subunit. No ligand binding is accompanied by a small leaking Ca^+ 
conductance. Free cytosolic proteins, nucleotides, and other substances can 
facilitate or impair IP3R gating by interacting with the IP3 binding-core. 



suppressor, and gate-keeper regions. Free endoplasmic-reticulum proteins 
and Ca^+ may also further modulate pore activity (not shown) via selectivity 
filters (small blue cylinders) located near pore helices (small rose cylinders). 
In the superposition state |i^), the IP3R samples all possible index values 
Tmax marking Dmax, returning an output x denoting a successful or 
unsuccessful search for shortest times or fastest rates m needed for 
store-released concentrations of free cytosolic Ca^+ to continuously diffuse 
and autocatalytically activate a nearest neighbor receptor. This superposition 
state may be regarded indefinitely stable in saturating IP3 concentrations. A 
phase shift by 0|p3r reversibly inactivates the receptor channel with high free 
cytosolic Ca^+ concentrations. Another subsequent phase shift reactivates 
the channel, confirming solution m has been found. See Equation 10 and 
relevant text for additional details. IP3R conformation representations 
adapted from Clark and Eisenstein (2013) with permission. 



be expected to sample and to then mark problem solutions 
through the unitary action of the IP3R Oracle: 

l^fXlO) - |l»/2^/^ -I'^'^'UXIO) - |l>)/2l''^ (10) 

where again \x} denotes the index-register qubit set to |0> (i.e., 
all if values, including value 0 for the inactivated IP3 -bound 
IP3R conformation) and(|0) - |l))/2^/^ is the superposed IP3R- 
Oracle qubit. The IP3R Oracle phase-shifts \x) only when 6(x) 
returns 1 as a result. Function 6{x), accordingly, may be thought 



of as the low- and high-affinity Ca + binding sites located at the 
cytosolic end of the IP3R emulating Grover's quantum algorithm 
(Figure 2). Recall that IP3R activity demonstrates a bell-shaped 
response profile to cytoplasmic Ca^+ levels. Ca^+ generally excites 
and blocks IP3RS in respective low (e.g., ~50 nM-1 [lM) and high 
(e.g., > 10 (xM) concentrations. Low Ca^+ concentrations capable 
of inducing sustained Ca^+ waves bind to high-affinity bind- 
ing sites, changing cation-pore permeability by repositioning the 
transmembrane gate and activating Ca^+ conductance. In con- 
trast, very low or very high Ca^+ levels keep the IP3R in a closed 
state. Ca^+ binding sites operating as 9 (x) return 1 upon ion-pore 
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opening and otherwise return 0 as a result. Coincident with 6 {x) 
action and all-or-none IP3R stimulation, the IP3R Oracle labels 
|x> by flipping the state from |1> to |0>. Isolating the Oracle's 
operation to a particular IP3R molecular region and function 
presents some conceptual difficulties. Arguments could be made 
for cytoplasmic and ER-lumen channel sensors which detect the 
presence of cations ions and IP3R interactions with nucleotides, 
proteins, and other substrate. Regardless, as noted above, once 
an IP3R becomes active, respective high quantal Ca^+ release 
exceeding 10 |xM near the opening of the same channel begins 
to cause reversible autoinhibition, presumably via low-affinity 
Ca^+ binding sites located proximal to IP3R cytoplasmic termi- 
nus. Inactivation may last longer than the open-channel period r . 
Such an effect is tantamount to the \x) phase shift produced by the 
Oracle of Grover's quantum algorithm (Figure 2). Furthermore, 
channel inactivation spans the entire Fmax index set, placing 
the IP3R into another superposition state |i/r>, as expected from 
application of the second Hadamard transformation. When IP3R 
inactivation is finally reversed to a fully excited Ca^+ -permeable 
state, the conditional phase shift, \x) —(—1)®'^' \x), has been 
performed for all basis states not equal to |0>. The last Hadamard 
transformation puts the register qubit into an equally weighted 
superposition for possible future Oracle summons. Since r can 
be accurately determined to be greater than the time of intersite 
Ca^^ diffusion for continuous waves (Strier et al., 2003; Solovey 
and Ponce-Dawson, 2010) and since the IP3R-Oracle phase shift 
of \x) only serves to emphasize timescale differences implicit in F 
of saltatory and continuous Ca^+ waves, the IP3R-mediated fire- 
diffuse-fire model simulating Grover's quantum algorithm will 
find target m indexed to with quadratic improvement in search 
efficiency. Together, steps 3 to 5 of Grover's quantum algorithm 
may be expressed as: 

H®"^(2|0>(0| -7)H®"-' =2|0>(0| -7, (11) 

where I is the identity matrix (cf. Nielsen and Chuang, 2000). The 
overall effectiveness of arriving at a solution involving a fast intr- 
acluster continuous Ca^+ wave predictably grants neurons better 
opportunities to initiate local and global response regulation for 
a variety of necessities. 

TESTABLE PREDICTIONS AND SIGNIFICANCE OF GROVER'S 
QUANTUM ALGORITHM FOR OTHER CICR MODELS 

A fundamental prediction of the above analytical treatment unad- 
dressed by standard classical interpretations is that, for an indi- 
vidual IP3R-based Grover's quantum algorithm to arrive at a 
solution m = Dq upon sensing/actuating F^ax with multiaffinity 
IP3R cytosolic Ca^+ binding sites and gating kinetics, quadratic 
and exponential speed-ups in respective algorithmic search time 
and wave velocity must coexist, with a maximum algorithmic 
search time of ©(Dmax/Dc)'''^ and wave velocity of vc = v^. 
If these constraints are not met, algorithm-processing capabili- 
ties will approach a classical algorithm taking 0(Dniax/£'s) time. 
After only quick inspection, the velocity equations for saltatory, 
vs = {D/d)g^^r, and continuous waves, vc = (D/r)'/^/^^F, 
may falsely imply to readers that a quantum-search result is 
impossible. Indeed, the arguments of Ponce-Dawson et al. (1999), 



for instance, enforce the idea that (global or intercluster) stable 
saltatory wave propagation is the fastest mode of transmission, 
at least for oocyte maturation. That conclusion heavily relies on 
the condition of rs < tc, where receptor-channel open dura- 
tion Tc might be one to two orders of magnitude larger than rj. 
This type of variation in receptor open time is not experimen- 
tally reported for IP3RS, which tend to be open for a fixed period 
between 15 and 40 ms depending on receptor subtype irrespec- 
tive of Ca^+ wave propagation mode (Foskett et al., 2007; Clark 
and Eisenstein, 2013). And, if t is instead calculated to be the 
rise time of Ca^+ concentration across a local grouping of chan- 
nels contributing to a release event, the divergence between rs 
and Tc can be expected to be no more than one magnitude for 
most neuronal Ca^+ waves. Taking these aspects into considera- 
tion and utilizing examples of experimentally realistic parameter 
values (e.g., Izu et al, 2001; Strier et al, 2003; Foskett et al., 
2007; Clark and Eisenstein, 2013) for d = 2 jxm between IP3R 
clusters, D = 190 (im^/s for continuous waves, D = 15 (jim^/s for 
saltatory waves, and r = 0.04 s for both continuous and salta- 
tory waves, it becomes apparent continuous waves (y8 = 1.9, v ~ 
(D/r)'/^ ~ 69|xm/s) can exceed saltatory wave (jS = 0.15, v ~ 
D/d ~ 7.5 |xm/s) velocity by greater than a power of 2. Despite 
being a simple example, these values underscore the plausibility of 
fast continuous Ca^+ waves displaying characteristics consistent 
with those predicted by the Grover's quantum-algorithm model. 
Importantly, the version of the fire-diffuse-fire model used in 
the present article assumes deterministic channel refractivity and 
instantaneous buffering for a without specification of Ca^+ store 
re-uptake or extracellular extrusion (cf. Ponce-Dawson et al., 
1999). Fast and slow Ca^+ buffering and Ca^^" extrusion and 
sequestration, such as that parameterized in stochastic mod- 
els (e.g., Coombes and Timofeeva, 2003; Coombes et al, 2004; 
Keener, 2006), will decelerate and even quench wave propa- 
gation, especially lower concentration saltatory waves. But the 
impact buffering has on wave-conduction modality and veloc- 
ity is further dependent on model selection. For instance, in 
the rapid-Ca^+ -buffering approximation with or without slow 
re-uptake (e.g., Strier et al., 2003), the fire-diffuse-fire model fur- 
ther becomes susceptible to error when describing saltatory-wave 
evolution. Continuous waves, however, largely overcome buffer- 
ing capacity, even more so when a massive wave-induction event 
occurs, such as calcium overload from extracellular sources or 
large ER mobilization. These sorts of effects predictably serve to 
increase F, effective diffusion, wave velocity, and thus fast stable 
Ca2+ wave search patterns at scales of receptor clusters or larger 
ER membrane surface areas supported by the Grover's quantum- 
algorithm model. In situations where waves are quenched or 
annihilated through Ca^+ buffers, re-uptake, or other factors 
(Keizer et al., 1998; Thul et al, 2007), the affected area will act as 
a wave guide to direct the (stochastic or deterministic) initiation 
and movements of subsequent waves, including planar, spiral, and 
oscillatory wave profiles capable of sending and storing distinct 
types of cellular information. With respect to these contexts (cf. 
Falcke, 2003a), future detailed numerical examination of param- 
eters is needed to identify the continuum limits of fire-diffuse-fire 
saltatory and continuous wave dynamics and their relation to 
expression of Grover's quantum algorithm at different scales of 
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CICR organization. As a special case of percolation Ca^+-release/- 
diffusion universality (cf. Timofeva and Coombes, 2004; Solovey 
and Ponce-Dawson, 2010), fire-diffuse-fire models conform to 
different spatiotemporal scales. The Grover's quantum algorithm 
model, in its conceptual form as a single IP3R operator, is best 
understood through local CICRs conducting at the scale of sin- 
gle receptor clusters, so that the model may reside in a quantum 
regime at both thermodynamic and informational degrees of free- 
dom. For local or intracluster wave propagation mode and veloc- 
ity, where d = 20 nm, the set of above (valid) values for global 
(i.e., intercluster) continuous and saltatory waves yield respec- 
tive local (intracluster) continuous waves of /J = 19000 with v ~ 
69 |xm/s and = 1500 with v ~ 19 |xm/s. The large values for /S 
indicate conditions supporting saltatory waves at a global scale 
do not necessarily ensure presence of saltatory waves for intra- 
cluster dimensions. Thus, although the algorithm searches for and 
finds optimal solution m, the boundary between local and global 
waves needs to be better resolved through numerical examina- 
tion of Fmax — the critical Grover's quantum-algorithm index of 
effective diffusion and solution m. 

Though a computationally tractable term containing many 
essential traits, employing the deterministic threshold-dependent 
Tmax unsatisfactorily weakens the explanatory power of the cur- 
rent model in regard to channel gating kinetics and quantum 
molecular action underlying formation and interrogation of 
superposition states, such as IV'^). Resolution of state bound- 
aries for the IP3R Grover's quantum-algorithm model therefore 
should be perfected by redefining Tmax and, naturally, constituent 
r to include buffering, re-uptake, and gating kinetics terms, 
as is done with more complex wave-evolution equations (e.g., 
Strier et al, 2003; Timofeva and Coombes, 2004; Thul et al., 
2007). Deterministic (DeYoung and Keizer, 1992) and probabilis- 
tic (Falcke, 2003b) mathematical models of single IP3R behavior 
offer richer accounts of the dynamic range expected for chan- 
nel activity, wave profiles, feedback control, and search patterns 
inherent in different channel reaction kinetics and structural con- 
figurations (Figure 2). Accordingly, Equation 3 may be restated, 
for benefit of future investigations, as: 

r = (Oi - <D2 - <t)3) / ([Ca2+]6o,, - [Ca2+]w,„) . (12) 

Variables Oi, <t>2, and <i>i quantify changes in point-source 
store-operated Ca^+ emission through respective release, re- 
uptake, and buffering kinetics relative to threshold activation 
kinetics, ([Ca^+lwt,. - [Ca^+Jwfj, with [Ca^+Jto,,., [Ca^+l^to = 
[Ca^^^lr — [Ca^^"]fc. The role of molecular dynamics in channel 
state is evident in each O (cf DeYoung and Keizer, 1992; Falcke, 
2003b); Oi = CER/c(i^MOpO + i^i)([Ca2+]ER - [Cfl^+Jc), $2 = 
fu[Ca2+]2,/[Ca2+]2 -1-^2, and $3 = Ji:B[Ca2+]c[B]c, where 
Cer/c is the ratio of Ca^+ in ER and cytoplasmic volumes, Fmo 
is maximum outward IP3R Ca^+ flux, is IP3R Ca^+ leak 
flux, fu is maximum ER Ca^+ uptake, po is the probability of 
IP3R-channel opening, A is the ER-uptake-activation constant, 
is the Ca^+ -buffer-binding constant, [Ca^+]ER and [Ca^+]cF 
are free ER and cytoplasmic Ca^+ concentrations, and [B]c is 
the free Ca^+ buffer concentration. While the complexity of 



these equations may be increased, addition of the open-channel- 
probability term, 

Po = [([Ca2+] [IP3] Kcai) / ([Ca2+] [IP3] + [IP3] Kca, 
+ Km.Kca, + [Ca^+j Km,) ([Ca^+j + JCc«a)]' - (13) 

where Kcai is the dissociation constant for IP3R channel inhibi- 
tion, -fT/psj and are IP3 dissociation constants, and Kca^ is the 
dissociation constant for IP3R channel activation, gives a much 
stronger framework for molecular considerations. For example, 
a consequence of the new T formulation is that the initialized 
index-register qubit becomes \x) = 1/4'/^ Z^MnLo 1 100)- Vector 
1 100) is the molecular configuration of composite IP3R binding 
subunits for each IP3R monomer (Mn) primed for channel acti- 
vation and Ca^+ conductance. The first, second, and third vector 
columns of |100) represent the respective (saturated) occupied 
IP3-biding subunit (i.e., 1), the unoccupied high-affinity Ca^+- 
binding subunit (i.e., 0), and the unoccupied low-affinity Ca^+- 
binding subunit (i.e., 0). When multiqubit-controUed function 
6{x) = 1 = 1 110), the high-affinity subunit binds Ca^+, making 
the IP3R permeable to cations. The IP3R Grover's operator, also 
acting as a multiqubit operator through coligand channel gat- 
ing kinetics, as earlier conjectured, then labels |x) by flipping the 
final two or all column values in conjunction with saturating- 
or subsaturating-IP3, Ca^+ -dependent channel inactivation. The 
operators's rate and probability of flipping between 0 and 1 
for each column and of finding target solution m is expectedly 
proportional to percent saturation of [Ca^+]c for high- and low- 
affinity Ca^+ -binding subunits and, therefore, proportional to 
index Fc, Dq, and vq- 

The IPsR-based Grover's quantum algorithm of CICR behav- 
ior and Ca^+ -mediated cellular response regulation thus imparts 
greater comprehensiveness than possible with traditional single 
receptor, intracluster, and intercluster models. For instance, above 
quantum depictions of the IP3R conformational state vector 
\x) uniquely permits study of computationally efficient subcel- 
lular superdense coding (Clark, 2010c), quantum learning and 
memory (Clark, 2010a,b,c,d, 2011, 2012a,b, 2013a; Liu et al, 
2013), quantum error diagnosis and correction (Clark, 2010c, 
2013a), and quantum encryption (Clark, 2010c, 2013a, in press). 
Recalling the discussion on the third Bell state, (|00)-|ll))/2'/^, 
the IP3R Grover's quantum-algorithm model may be extended to 
quantum coupling between two nearest neighbor IP3 -saturated 
activation-primed receptor channels Ri and _R2- For Ri and _R2> 
each in superposition state \x) = 1 /4'/^ X!Mn=o 1 100), the entan- 
gled bell bases may be deemed bidirectional coupling factors, 
such as Ca^+ sensitivity, shared between receptors and imposed 
by level of cytoplasmic modulator (e.g., ATP) saturation on 
coligand-dependent IP3R allostery (Foskett et al, 2007; Clark 
and Eisenstein, 2013), where Ri and -R2 each possess one unique 
highly correlated complimentary Bell-state qubit. _Ri, depending 
on the action of multiqubit-operators 0 and Grover's iteration, 
may send classical information intrinsic to F to R2 in the form 
of two-bit strings, such as 00 indicating neither sets of Ca^+ 
binding sites are occupied, 10 indicating only the set of high- 
affinity Ca^+ binding sites are occupied, or 0 1 indicating only the 
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set of low-affinity Ca^+ binding sites are occupied. By employ- 
ing the decoding gate matching Ri's classical signal, Ri can 
decode the superdense signal sent from Ri with less computa- 
tional expenditure and superior acquired content than possible 
with classical information processing (cf Nielsen and Chuang, 
2000)\ Not only does this interpretation of channel conforma- 
tional state and channel-channel communication agree with the 
IP3R Grover's quantum algorithm and CICR cytophysiology, it 
also establishes conditions for, among other phenomena, the 
expression of quantum-molecular memory storage and retrieval 
at levels of single IP3RS and IP3R clusters. Since only three of 
four IP3R monomers must be activated for channel opening, 
the revised definition of |x> implies, in a manner like Ventura 
and Martinez (1999), single IP3RS must maintain or remem- 
ber Mn\/{Mn — 1)! = 4 superposed combinations of: (1) initial 
activation-primed conformational patterns 1 100) for about 19 ms 
or less to finish -labeling before IP3 dissociates from its IP3R 
binding subunit, (2) -transformed conformational patterns 
1 110) for about 600 ms or less to finish the Oracle call before 
Ca^+ dissociates from its high-affinity IP3R binding subunits, and 
(3) Oracle-transformed conformational patterns |101> or |001> 
for about 5 s or less to finish the algorithm's final phase shift 
before Ca^+ dissociates from its low-affinity IP3R binding sub- 
units (DeYoung and Keizer, 1992). These superposition states, 
embedded in |i/f), have spatiotemporal estimates well within 
quantum-decoherence bounds calculated for macromolecules 
located in live cells constrained by physiologically salient envi- 
ronments (cf Gutin et al, 1996; McFadden and Al-Khalili, 1999; 
Cieplak and Hoang, 2003; Davies, 2004) and are experimentally 
testable by molecular dynamics and CICR simulation as well as 
bioassays involving wildtype and selectively mutated IP3R iso- 
forms reconstituted in planar lipid bUayers. Although coverage of 
larger memory structures, such as that formed by multiply cou- 
pled IP3R intracluster patterns, is beyond the scope of this article, 
the same concepts presented for individual IP3RS are scalable to 
intracluster dimensions and organization. 

RELEVANCE OF GROVER'S QUANTUM ALGORITHM FOR 
HEALTHY AND DISEASED NEURONS 

In preceding sections, a fire-diffuse-fire model capable of explain- 
ing intracluster activity of individual IP3RS was identified and 
analytically parameterized as a candidate mechanism for a natural 
neuronal form of Grover's quantum algorithm. Model accu- 
racy fundamentally depends upon the sensitivity of IP3RS to 
physiological parameters characterizing Ca^^ -channel molecular 



^In this example of Bell state k)10 = (|00) -|- \n))/2^l^, the first |00) and 
second |11) qubits respectively denote equilibrium and high Ca^+ sensitivity 
induced by saturating and supersaturating free cytoplasmic ATP. JJj, wrhich 
possesses the first qubit of coupling factor k)io , applies 0(x) = \ = \x) = 
1 110), interacts \x) with its half of the coupling factor, and transmits clas- 
sical bit-string 10 via T. R2, which possesses the second qubit of coupling 
factor \K:)ig, will next employ its half of the coupling factor and associ- 
ated quantum-flip gate, S{x) = 0 — J- 1 = |x) = 1 1 10), to recover information 
about the conformational state |x> of Ri and initiate its own Grover's opera- 
tor. This computational scenario, resilient to error, is equivalent to three-agent 
quantum teleportation transpiring over classical communication channels (cf 
Nielsen and Chuang, 2000). 



Structure and function as well as scalable quantum-level gains in 
classical Ca^+ diffusion rates, Ca^+ wave propagation, and appro- 
priate fast cellular response regulation. Computations made by 
the IP3R algorithm infer search target solutions for fast classi- 
cal Ca^+ diffusion rates via interrogation of index variable Fmax 
associated with detectable free cytosolic Ca^+ concentrations and 
properties of IP3R channel conductance. Collectively, these fea- 
tures of the model allow the algorithm to operate within quantum 
computational and thermodynamic regimes without concern of 
incurring statistical mechanics measurement problems, such as 
decoherence of processed superposed eigenstates. The model 
nevertheless yields only a computational first-approximation of 
Grover's quantum algorithm and needs future refinement by 
applying sophisticated: ( 1 ) relativistic quantum physicochemistry 
theory to aptly match IP3R protein structure and function with 
the workspaces and operators of Grover's quantum algorithm, 
and (2) fire-diffuse-fire or lattice-percolation mathematical treat- 
ments of intracluster IP3R activity and CICR dynamics to fully 
address aspects of reaction-diffusion stochasticity and cytosolic 
Ca^+ buffering (cf Izu et al., 2001; Strier et al, 2003; Guisoni 
and de Oliveira, 2006; Solovey and Ponce-Dawson, 2010). With 
respect to the latter topic, the present model becomes partic- 
ularly relevant during actuation of intracellular compartmental 
Ca^+ loading from interstitial and intracellular cation sources. 
Even as a preliminary construct, the model implies contexts coin- 
cident with moderate to massive fluxes of Ca^^ through cation- 
permeable integral cell membrane pores and gated channels, such 
as during synaptic plasticity (Malenka and Bear, 2004), microbial 
pathogen attack (Clark, 2013b; Clark and Eisenstein, 2013; Clark 
et al, 2013), pathological oxidative stress (Benedicte et al, 2012; 
Clark, 2012c), and neurological disease and aging (Verkhratsky, 
2005; Bezprozvanny and Mattson, 2008; Stutzmann and Mattson, 
2011), will assist in driving neurons to accelerate response regula- 
tion to quantum-level efficiency through induction of stable local 
and possibly subsequent global continuous Ca^+ waves. From a 
physiological perspective, dramatic increases in Ca^+ wave veloc- 
ity and signal transduction at either intracluster or intercluster 
physical dimensions are impressive and attainable for a all sorts 
of differentiated eukaryotic cells (cf Izu et al., 2001), requir- 
ing a maximum algorithmic search time of ©(Dmax/Dc)'^^ and 
wave velocity of vc = to be realized. However, whether it is 
activated by stochastic blips of high Ca^+ -conductance or deter- 
ministic cellular Ca^+ loading, a neuronal version of Grover's 
quantum algorithm, just like that proposed for other eukary- 
otic cells, figures to promote advantageous subcellular superdense 
coding (Clark, 2010c), quantum learning and memory (Clark, 
2010a,b,c,d, 201 1, 2012a,b, 2013a; Liu et al, 2013), quantum error 
diagnosis and correction (Clark, 2010c, 2013a), and quantum 
encryption (Clark, 2010c, 2013a, in press). 

Equally significant, the quantum computational value of a 
receptor-scale Grover's quantum algorithm also can be expected 
to contribute to surprising classical information processing over 
much longer intracellular distances and times common to global, 
multicompartmental Ca^+ signaling. This expectation, unac- 
counted for by standard (stochastic or deterministic) CICR mod- 
els, agrees with expression of three-agent quantum teleportation 
over communication channels transmitting classical bits via the 
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information content of F, a circumstance, as described in the 
above section, that exposes superdense coding and quantum 
molecular memory for IP3R systems. Limited evidence suggests 
digital representation of cellular processing is encoded, transmit- 
ted, and stored by free intracellular Ca^+ (Plieth, 2005), CaMKII 
holoenzyme (Hameroff et al., 2010), and other Ca^^-related 
substrate broadly distributed across affector-effector systems. 
Individual molecules conveying bitwise information may form 
higher-order bytes at large concentrations and/or with molecular 
complexes, such as that reported for six-domain CaMKII encod- 
ing of microtubule lattices (Hameroff et al., 2010). Expression 
of an IPsR-based Grover's quantum algorithm, which may serve 
as a quantum amplifier and router, supports an estimated forty- 
fold boost in classical information processing by networked Ca^+ 
release sites through buffer-dependent superadditive Ca^^ wave 
densities and velocities (c£ Izu et al, 2001). Such effects presum- 
ably enhance operational traits of, for example, signal coincidence 
detection and integration, bidirectional synaptic plasticity, gene 
expression, immunodefenses, growth and tropisms, protein mod- 
ification and transport, cytoskeletal polymerization, endosome 
formation and other cell functions by rapidly selecting, order- 
ing, and/or counting optional local response regulation choices. 
The impact on neuronal cyctoskeleton operation alone attests 
to favorable cascading effects governing geometry of dendritic 
spines and synaptic cleft widths, intracellular molecular and 
vesicular trafficking, membrane repair, synaptogenesis, neurite 
growth, and efficacious synaptic transmission (e.g., Malenka and 
Bear, 2004; Verkhratsky, 2005; Bezprozvanny and Mattson, 2008; 
Craddock et al, 2010; Priel et al, 2010; Dent et al, 2011). 
Using archetypal glutamatergic neurons (e.g., Verkhratsky, 2005; 
Hagenston and Bading, 2011), which are exquisitely sensitive to 
fluctuations in intracellualar Ca^+ homeostasis, one can readily 
extrapolate how the IPsR-based Grover's quantum algorithm fig- 
ures to help execute ER-dependent signal amplification and inte- 
gration in healthy cell states of all neurons (Clark, 2012b). Unlike 
spontaneous punctate store-operated Ca-^+ emissions observed 
for every eukaryotic cell or environmentally triggered massive 
store-operated Ca^+ overload observed for certain cell types 
(e.g., oocytes), large temporary elevation of Ca^+ microdomain 
concentrations following post-synaptic Ca^+ entry through acti- 
vated NMDARs initiates widespread and differential response 
regulation in dendritic, somal, and axonal compartments. The 
amount and spread of inward Ca^+ current is too small to 
effect most transduction processes. Instead, NMDAR-mediated 
Ca^"*" entry stimulates secondary Ca^+ release from intracellu- 
lar stores. Information processed by NMDAR-dependent Ca^+ 
cascades is augmented by glutamate diffusion to extrasynaptic 
GPCRs responsible for intracellular IP3 generation. Compared 
to RyRs, IP3RS tend to have greatest density in ER membranes 
located in the soma and dendritic shafts of neurons, where col- 
igand IP3 may exert maximum influence (c£ Stutzmann and 
Mattson, 2011). Thus, an IP3R Grover's quantum algorithm 
likely guides heterosynaptic activity as well as kinase-induced 
(e.g., CaMKII and IV) gene transcription and protein synthe- 
sis accompanying longer-term structural plasticity, chiefly that 
of LTP (cf. Fitzjohn and CoUingridge, 2002), rather than vesic- 
ular events particular to shorter-term presynaptic paired-pulse 



and post-tetanic facilitation. The algorithm ensures synaptic plas- 
ticity maintenance and nuclear response regulation by selecting 
the best target solutions for fast classical Ca^+ diffusion rates, 
stable continuous Ca^+ wave modes, and, consequently, activa- 
tion of kinase and other messenger pathways via detectable free 
cytosolic Ca^+ concentrations and properties of IP3R channel 
conductance. 

However, the appearance of a quantum-efficient search algo- 
rithm in the function of neurons or, for that matter, any other 
cell type need not necessarily guarantee evolutionary and/or eco- 
logical benefit for the cell that implements the algorithm and the 
host to which the cell belongs. Indeed, viral, bacterial, fungal, and 
protozoal infectious agents, including, among other pathogens. 
Human Immunodeficiency Virus type 1 and Trypanosoma par- 
asites, that selfishly coopt neuron intracellular Ca^+ systems 
can use the algorithm to optimize timing and effectiveness 
of infection stages against barriers to invasion, pathogenesis, 
proliferation, and release. While most pathogens deploy well- 
timed Ca^+ -dependent trophic and deleterious strategies, such as 
genomically encoded proteins and lipopolysacchrides, to exploit 
host-cell physiology, infected neurons are unusually susceptible 
to metabolic distress, apoptosis, and additional harmful effects 
leading to host cognitive impairments (Clark, 2013b; Clark and 
Eisenstein, 2013; Clark et al., 2013). Similarly, an uninfected, 
but diseased or aging, neuron can diminish its own cell perfor- 
mance by speeding-up the selection and execution of cellular 
response regulation incompatible with cell or host survival. Good 
examples involve neurodegenerative Parkinson's and Alzheimer's 
diseases, where ER-mediated deficits play major roles in disease 
severity and progression (cf. Stutzmann and Mattson, 2011). 
In the case of Alzheimer's disease, persistent upregulation of 
intracellular Ca^+ concentrations are known to commence and 
accelerate synaptic loss and amyloid plaque disposition, manifest 
pathocytological characteristics of the disease. While the excito- 
toxic nature of Ca^+ dysregulation in Alzheimer's disease is well 
understood, the exact causes are not (Stutzmann and Mattson, 
201 1; Popugaeva and Bezprozvanny, 2013). Recent and somewhat 
controversial evidence indicates that IP3R activation by prese- 
nilins, protein products of autosomal inherited mutated genes 
PSl and PS2 linked to early onset Alzheimer's disease, may be a 
major contributing factor. The corresponding heightened IP3R- 
mediated Ca^+ mobilization means implementation of IP3R 
Grover's quantum algorithm drives diseased, injured, and/or 
aging neurons toward faster catastrophic failure than would be 
otherwise possible with classical response regulation. Hence, the 
delicate balance between beneficial and detrimental uses of a neu- 
ronal version of Grover's quantum algorithm demands thorough 
theoretical and empirical scrutiny of cellular conditions gov- 
erning the algorithm's application in both healthy and diseased 
states. 
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